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Abstract 



The unification of quantum field theory and general relativity is a 
fundamental goal of modern physics. In many cases, theoretical efforts 
to achieve this goal introduce auxiliary gravitational fields, ones in ad- 
dition to the familiar symmetric second-rank tensor potential of general 
relativity, and lead to nonmetric theories because of direct couplings be- 
tween these auxiliary fields and matter. Here, we consider an example of 
a metric-affine gauge theory of gravity in which torsion couples nonmin- 
imally to the electromagnetic field. This coupling causes a phase differ- 
ence to accumulate between different polarization states of light as they 
propagate through the metric-affine gravitational field. Solar spectropo- 
larimetric observations are reported and used to set strong constraints on 
the relevant coupling constant k: k 2 < (2.5 km) 2 . 

1 Introduction 

The quest for a complete and self-consistent unification of quantum field theory 
with the theory of general relativity remains one of the most important unsolved 
problems in modern theoretical physics. Its solution promises an extension 
of the standard model of particle physics that encompasses all fundamental 
interactions and accounts for everything from the evolution of the early universe 
to black hole physics. This possibility continues to inspire work on such theories 
even after decades of effort. Clearly, one can expect candidate theories to predict 
new physics, that is, phenomena beyond the scope of general relativity and the 
standard model. The birefringence we study here is an example. 

Quantum gravity research is unusual in that it is driven by the need to over- 
come extraordinary conceptual and mathematical difficulties in formulating a 
consistent theory rather than by an accumulation of experimental evidence that 
could inform ones choices among theoretical alternatives. Exploring such alter- 
natives during the past several decades has produced many theories of gravity 
and, thus, a need for an overarching theoretical framework within which these 
alternatives can be compared and classified. The Dicke framework, defined in 
Appendix 4 of Dicke's 1964 Les Houches lectures [I], can be seen as the basis for 
a number of formalisms encompassing Lagrangian-based local field theories of 
gravity. It demands that gravity be associated with one or more fields of tenso- 
rial character. Consequently, alternatives to general relativity generally feature 
one or more gravitational fields in addition to the usual second-rank symmetric 
tensor potential. The distinction between metric and nonmetric alternatives 
often turns on the manner in which matter couples to such additional fields. 
The purely geometrical character of general relativity and other metric theories 
is a consequence of the familiar minimal coupling of matter fields to a single 
symmetric second-rank tensor gravitational field. Nonmetric theories deviate 
from this pattern, typically including direct couplings between matter and aux- 
iliary gravitational fields. The Dicke framework encompasses nonmetric as well 
as metric theories of gravity and, therefore, general relativity as a particular 
example for a metric theory. In this paper we focus on observable consequences 
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of a nonminimal coupling between a gravitational torsion field and the electro- 
magnetic field. 

An appropriate framework for the analysis of electrodynamics in a back- 
ground gravitational field is given by the xg-formalism, invented by W.T. Ni 
in 1977 2 . The \ of its name refers to a tensor density which provides a 
phenomenological representation of gravitational fields. The structure of the 
Xg-formalism is in agreement with the basic assumptions and constraints of 
the Dicke framework. Demanding electromagnetic gauge invariance and linear- 
ity of the electromagnetic field equations, the most general Lagrangian density 
governing source-free electromagnetic field dynamics is 

£ng = -^x aplS F a ^ 5 . (i) 

The independent components of the tensor density ^"^T" 5 comprise 21 phe- 
nomenological gravitational potentials that allow one to represent gravitational 
fields in a very broad class of nonmetric theories. In 1984 Ni [3j noted that the- 
ories encompassed by this formalism can predict birefringence and used pulsar 
polarization observations to constrain this possibility. Perhaps because no spe- 
cific theory predicting birefringence was recognized at that time, the significance 
of Ni's result was initially overlooked. 

This situation changed when J.W. Moffat published a revised version of his 
non-symmetric gravitation theory (NGT) [5] in 1990. One year later Gabriel 
et al. showed that NGT predicts a polarization-dependent speed of light as a 
consequence of a violation of the Einstein equivalence principle and imposed 
the first sharp constraints on the magnitude of this birefringence |H1 E] ( m 
this context see also Solanki & Haugan 1996 @J and Solanki et al. 1999 [H])- 
In contrast to these earlier investigations, the present paper is not concerned 
with NGT. It focuses on another alternative to general relativity drawn from 
the class of metric-affine theories of gravity (MAG) m part, because an 
array of technical difficulties have challenged the viability of NGT We 
show that birefringence is a generic feature of MAG in the case of nonminimal 
coupling between electromagnetism and gravity, specifically torsion, and use 
solar data to set sharp constraints on a relevant coupling constant. We hope 
that this example will prompt further work on gravity-induced birefringence 
predicted by metric-affine theories and other alternatives to general relativity. 
In that regard, see ^] El • 

It is important to note that the significance of such constraints goes far be- 
yond testing versions of NGT or MAG. These are merely concrete examples 
of theories that predict gravity-induced birefringence, a phenomenon shown by 
Haugan and Kauffmann 10 to be predicted by a broad class of the nonmetric 
theories encompassed by the x<?-formalism. Here we make use of the formalism 
they invented to compute the effect of gravity-induced birefringence using the 
X9 representation of any gravitational field. We emphasize that observations 
constraining the strength of such birefringence complement more familiar tests 
of the Einstein equivalence principle like the Eotvos, gravitational redshift and 
Hughes-Drever experiments, which currently do not bound our special nonmin- 
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imal coupling. The reason is that we focus on a quite novel coupling between 
electromagnetic fields and torsion which was not taken into account in the anal- 
ysis of these classical experiments. However, since our new coupling also violates 
local Lorentz invariance as one can see from Eq. (J2J, it is indeed thinkable that 
Eotvos, Hughes-Drever or redshift experiments could also constrain this new 
possibility. 

We also report on new solar polarization data and the constraints on such 
gravity-induced birefringence obtained therewith. In section III. A we describe 
the new profile difference technique we use to search for evidence of birefrin- 
gence. In section III.B we review the Stokes asymmetry technique proposed 
previously by Solanki and Haugan (1996) 4 . In section IV we describe the 
solar data we analyze and the observations that produced them. The analy- 
sis of these data and the constraints we infer on gravity-induced birefringence 
are discussed in section V. We quote constraints on the Sun's NGT charge, Iq, 
purely as a figure of merit that can be compared to prior constraints, but focus 
on constraints on a metric-affine parameter, fc 2 , defined in the next section. 

2 Gravity-induced birefringence in metric-affine 
theories 

An empirically adequate account of gravitation is given today by general rel- 
ativity which predicts vanishing torsion and vanishing nonmetricity (covariant 
derivative of the metric). This account is, so far, in complete agreement with 
observational results. However, general relativity is a classical theory. The de- 
sire for a quantum mechanical account of gravitation requires a more general 
framework. 

One possibility is to formulate gravity as a gauge theory of an underlying 
local spacetime group. Metric-affine gravity (MAG) is based on the assumption 
that affine transformations are the gauge group. It is the most general canonical 
gauge theory of gravity [T3] . 

Metric-affine theories use a second-rank symmetric tensor field, a co-frame 
field and a connection one- form field to represent gravitational potentials. Al- 
though the symmetric tensor is referred to as the metric, metric-affine theories 
are generally nonmetric. Other gravitational potentials, specifically the torsion 
and nonmetricity fields extracted from the connection, couple directly to mat- 
ter. Nonmetric couplings to the electromagnetic field can lead to gravity-induced 
birefringence. 

In this section we will carefully examine the possible coupling of the elec- 
tromagnetic field to gravity within the framework of MAG. Under the assump- 
tion of electric charge and magnetic flux conservation Puntigam et al. [T?)] 
showed that the conventional formulation of Maxwell-Lorentz electrodynamics 
in a metric-affine gravitational field predicts that the propagation of light is not 
influenced by the presence of post-Ricmannian geometric fields like torsion T a 
or nonmetricity Q a p. However, "admissible" nonminimal coupling possibilities 
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that lead to birefringence do exist since one can modify the Maxwell-Lorentz 
spacetime relation H = X *F. Couplings of this kind respect gauge invariance 
and, as a consequence, electric charge as well as magnetic flux conservation. 
In order to make quantitative predictions about electromagnetic field dynamics 
in a metric-affine background field using this nonminimal coupling scheme, a 
specific additional Lagrangian density is needed. We consider 

L EM = k 2 * (T a A F) T a A F, (2) 

consistent with the suppositions made above. Here, A: is a coupling constant 
with units of length, * denotes the Hodge dual, T denotes the torsion and F 
the electromagnetic field. Our intention is to set strong limits on k 2 and, so, to 
decide about the physical relevance of the gravity-induced birefringence gener- 
ated by the coupling J5J. Currently there are no expectations from the theory 
side of what a reasonable value for k 2 might be. The additional Lagrangian J5J 
can be written as 

5C EM = 5x aM5 F a pF lS , (3) 

so that we can make use of the general formalism developed by Haugan & 
Kauffmann in ^U] to infer the consequences of the nonminimal coupling J5J 
on the propagation of light through a metric-affine field. This means that our 
objective is to compute the fractional difference in the propagation speed of 
linear polarization states 

^ = \V(A-Cy + 4B2 . (4) 

that are singled out by a solar torsion field. The coefficients A, B and C depend 
on the location in space-time and on the direction in which the wave propagates. 

Since we are going to search for birefringence in the essentially static, spher- 
ically symmetric gravitational field of the Sun, we are interested in static and 
spherically symmetric solutions of the metric-affine field equations for torsion. 
One such solution was given by Tresguerres in 1995 ^[^H- It can be split into 
nonmetricity dependent and independent parts. The latter, which is assumed 
to couple to the electromagnetic field via (0, is given by 

A 9 a , (5) 

with torsion mass m and with kg = 1 for a — 0, 1 or fco = — 1/2 for a = 2, 3, 
respectively. In this equation we have dropped the dilatation charge No from 
Tresguerres' original solution because it vanishes if the nonmetricity field does. 
Also the small observed value of the cosmological constant A means that effects 
of its term can be neglected on galactic and smaller scales. 

Equation JSJ is consistent with the most general static, spherically symmetric 
form for a torsion field ^7] 

T° = a(r) 9° A 9 1 + a{r) 9 2 A 3 , (6) 
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T 1 = (3{r) 9° Ad 1 + f3(r) 9 2 A 9 3 , (7) 
T 2 = 7(1) 9° A 9 2 + 7(2) 9° A 9 3 

+ 7(3) 9 1 A 9 2 + 7(4) 1 A 9 3 , (8) 
T 3 = 7(1) 9° /\ 9 3 ~ 7 (2) 6»° A 9 2 

+ 7(3) 1 A 3 — 7(4) 9 1 A 2 . (9) 



The solution JSJ is a special case having <5(r) = j3(r) = 7( 2 ) = 7(4) = 0. 

Plugging this general representation of a spherically symmetric torsion field 
into the Lagrangian density yields 



SC EM = k 2 {(a 2 - P 2 )B( - (7^ + ^)[Bi + B 2 } 

- (7(3)+7( 2 4))[^2 2 +^3] (10) 

+ 2 (7(i)7(4) - 7(2)7(3) )[B 2 E 2 + B 3 E 3 ] 
+ 2 (7(i)7(3) + 7(2)7(4) )[B 3 E 2 - B 2 E 3 }} , 

where E and B refer in the usual way to the electric and magnetic components 
of F in Tresguerres' coordinates. Note that in his notation the three-vector 
index 1 refers to the radial direction. In terms of the components of the £, £ 
and 7 tensors of reference |10j . the corresponding Lagrangian expression is 

5C-EM = C11B 2 + ( 22 B 2 + (33-63 _ £22-62 - C33E 2 (11) 
+ 2{j 22 E 2 B 2 + j 33 E 3 B 3 + -f 32 E 3 B 2 + r y 23 E 2 B 3 }, 

from which the expressions for the nonzero components of the £, C and 7 ten- 
sors in terms of Tresguerres' functions can be read for the nonminimal torsion 
coupling considered here. Haugan & Kauffmann |1(J| have shown that the ex- 
pressions A — C and B from (£Q) can now be expressed in terms of the spherical 
components of these tensors £, Q and 7 

a - e = * + ei) + ^a - 7^) + (d? + a)) (12) 



v / 6 V 



and 



b = -±(i{&> - ao + 2(a + Tiio + - ao) . (13) 

Only / = 2 components appear. The m' notation indicates that these spher- 
ical components correspond to Cartesian components in a quasi-Lorentzian 
(t, x' ,y' , z') coordinate system oriented so that the radiation of interest prop- 
agates in the z' direction. In the following, we exploit the fact that these 
spherical components can be expressed in terms of components in another quasi- 
Lorentzian (i, x, y, z) coordinate system via the familiar transformation law 
e.g. 

$=Z$ TO fo^)&° - (14) 
where <fi, 9 and ip are the Euler angles specifying the rotation from (t, x, y, z) to 
(t, x' , y' , z') and the rotation matrix T> is given in terms of spherical harmonics. 
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Since the only direction that a spherical field can single out is the radial 
one, it is now useful to introduce, at each point along a light ray of interest, a 
local quasi-Lorentzian (i, x, y, z) coordinate system oriented so that the z axis 
is radial and the x axis lies in the ray's plane. This is convenient because the 
spherical tensors introduced above are simple in such local coordinate systems. 
Specifically, , Q? and 7™^ are nonzero only for m = 0. 

At each point along the ray, the local (t, x, y, z) coordinate system is rotated 

about the y axis through an angle 9 to obtain a local (t, x' , y\ z') system so that 

now the ray runs in the z' direction. The local value of Sc/c in Q is expressed 

(2) (2) 

in terms of A — C and B which are, in turn, expressed in terms of £±2'' C±2' 
and 7^]/ according to l|12fl and Since the Euler angles of the rotation 

from (t, x, y, z) to (t, x', y', z') are 9 and <j> = %j) = 0, the transformation law l|14fl 
simplifies to 

^l=Bm 2 e^ , (15) 

with the same relationship between q5, and Co 2 '' and between 7^, and . 
Here, 9 denotes the angle between the light ray's propagation direction and 
the radial direction. Exploiting these transformations, Haugan & Kauffmann 
conclude from (|12|l and 113|) that 

A-C = -j=($ ) +tf ) )sin 2 fl (16) 

and 

e = --L( 7 ( 2 ) sinfl) 2 . (17) 

Substitution into yields an expression for the fractional difference in the 
propagation speeds of linear polarization states for the case of a static spherically 
symmetric torsion field 



H } ) ■ (is) 



To express this in terms of Tresguerres' parameters we simply use the relation- 
ship between the spherical and Cartesian components of £, ( and 7 and between 
those Cartesian components and Tresguerres' parameters. 



So 


= e 1 


+ ^ 22 + f 3 = 


2fc 2 (7(3) +7(4)) 


A2) 

So 


= e 1 




) = -fc 2 (7( 2 3) +7 ( 2 4)) 


AO) 
^0 


= c 11 


+ (22 + C 33 = 


fcV-/? 2 ) -2(7(1) +7(2))] 


A*) 
^0 


= c 11 


_1 (C 22 +C 33 


)=fc 2 [( a 2 -/3 2 ) + (7 ( 2 1) +7( 2 2 ) )] 


(0) 

To 


= 7 11 


+ 7 22 + 7 33 = 


2fc2 (7(1)7(3) -7(2)7(4)) 


(2) 

To 


= 7 U 


_l( 7 22 + 7 33 


') = - fc2 (7(i)7(4) - 7(2)7(3)) ■ 
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For the particular torsion field (JSJ we have 7(2) = 7(4) = and, therefore, 
B = 0. Using JSJ) one gets 



> (2) -]/."";/' k 2 f 1 2m 



which yields 



so that we have 



4 

#, + tf, = _3*m _ (21) 



^-C = --^_sin 2 . (22) 

Therefore, the fractional difference between the velocities of the two polarization 
states is given by 

8c ,-k 2 m . 9 „ .„„. 

— = — \/6— 5- sm 2 9 , (23) 
c r } 

where 6 again denotes the angle relative to the outward radial direction in 
which the light is propagating. For the total phase shift A<I> which accumulates 
between the source and the observer one now has to calculate 

— dt = -V6ujk 2 m — —dt . (24) 



c 

Recapitulating the analysis performed in |3IHj , we find 

[2 2^(^ + 2)^-1) 
V 3 \R 2 m + 1 

where /i denotes the cosine of the light source's heliocentric angle, A is the light's 
wavelength and R is the Sun's radius. 



3 Technique 

We follow two strategies to test for gravitational birefringence. One of these 
was outlined by Solanki & Haugan (1996) 4 , but could not be applied due to a 
lack of appropiate data. It is summarized and its implementation is described in 
Sec. III.B. The other technique is new and is described in Sec. III. A. In order to 
compare the effectiveness of these techniques with results of previous attempts 
to set limits on gravity-induced birefringence, which dealt with NGT, we also 
briefly consider NGT here in addition to metric-afhne theories of gravity. 
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3.1 Profile Difference Technique 

This technique relies on the fact that A<& is expected to be a strong function of 
it, which is confirmed in the two concrete cases of NGT (see reference [§]) and 
metric-affine theories (see Sect. II). A$ is the phase shift which accumulates 
as light propagates from a point on the solar surface to the observer between 
net circular polarization, described by the Stokes parameter V, and net linear 
polarization aligned at 45° to the nearest part of the solar limb, generally as- 
cribed to Stokes U. Formulae for A<i> as predicted by metric-affine theories are 
given in Sect. 2. For Moffat's NGT (Moffat 1979 6 ) a corresponding expres- 
sion has been published by Gabriel et al. (1991)[7| This means that for 
any sufficiently large NGT charge £q or equivalent metric-affine parameter k a 
mixture of Stokes V and U profiles will be observed from sources of polarization 
distributed across the solar disc (i.e. as a function of fi), irrespective of the 
exact polarization state of the emitted radiation (i.e. which mixture of V and 
U is produced at the solar surface). Let the subscripts 'src' and 'obs' signify 
the Stokes profiles as created at the source and as observed, respectively. Then, 
irrespective of the value of < | V^ rc | > — < |J7 src | > for sufficiently large £q or 
k, < \V Q b s \ > — < |C/ bsl > will t en d to zero. The averaging is over all fi values 
and the total number of profiles is assumed to be large. 

This effect is illustrated in Fig.^ In Figs.^i & b iV^^sl — l^obsl is plotted 
vs. k and /z for the extreme cases \V SVC \ — (Fig. |TJi) and \U src \ = (Fig.^ 3 ). 
Other combinations of |Vst C (m)I an d \U src (ii)\ give qualitatively similar results. 
The |V^bsl — l^obsl surface oscillates ever more rapidly with increasing k and 
with decreasing /i. Lines of equal | V^bs I — \U bs\ are strongly curved in the k — fx 
plane. These two points combine to lead to decreasing < | V^bs I > — < |tobsl > 
with increasing k. This is shown in Figs. ^ & d for the cases illustrated in 
Figs, nji and b, respectively. As expected, the < | Vobs I > — < l^obsl > vs - k 
curves exhibit a rapidly damped oscillation around zero. This effect can be used 
to set upper limits on gravitational birefringence if the observations exhibit a 
< | Vobs I > — < | Cobs I > that differs significantly from zero. As we show later 
in this paper, this is indeed the case. 

3.2 Stokes asymmetry technique 

The strategy proposed by Solanki & Haugan (1996) .4: makes use of the symme- 
try properties of the Stokes profiles produced by the Zeeman splitting of atomic 
spectral lines. In the absence of radiative transfer effects in a strongly dynamic 
medium net circular polarization, Stokes V , is antisymmetric in wavelength and 
net linear polarization aligned at 45° to the solar limb, Stokes U , is symmetric. 
Birefringence, predicted by metric-affine gravity theories (and NGT) changes the 
phase between orthogonal linear polarizations and thus partly converts Stokes 
V into Stokes U and vice versa. However, U produced from V by gravitational 
birefringence still has the symmetry of V and can thus be distinguished from 
the Zeeman signal. 

Let subscripts 'a' and V signify the antisymmetric and symmetric parts of 
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the Stokes profiles, respectively. Then, 



17. 



a,src 



s,src 



s,src 



K.obs sin A$ + U^ ohs cos A$ 
y S;0bs sin Ai> + C/ Si0bs cos A$ 
^s.obs cos A$ + U sobs sin A$ 



(26) 



(27) 



la 



a,src 



V; i0bs cos A$ + J7 aj0bs sin A$ 



Thus for observed symmetric and antisymmetric fractions of U and V Eqs. (3.1) 
and lf2*?|) predict the ratios Ua,/U s and 14/14 at the solar source. 

If the solar atmosphere were static these ratios would vanish (L4,src = 
Va,src = 0), so that any observed 14 or V a would be due to either A$ or 
noise: £4 j0 bs = 14 jS rc sin A<I>, V s obs = /7 s ,src sin A$. The solar atmosphere 
is not static, however, and consequently the Stokes profiles do not fulfill the 
symmetry properties expected from the Zeeman effect even for rays coming 
from solar disc centre (/i = 1), which are unaffected by gravitational birefrin- 
gence. This asymmetry has been extensively studied, in particular for Stokes 
V, which most prominently exhibits it |2(JI 1211 f?2l I24j . Although most profiles 
have 14/14 < 0.2, a few percent of V profiles exhibit V B /V a values close to unity, 
even at \i = 1. Such profiles occur in different types of solar regions, e.g. the 
quiet Sun [221 > active region neutral lines [23 and sunspots j2E| ■ The magnitude 
of Vs/Vg, decreases rapidly with increasing Va and profiles with V s /Va, > 1 are all 
very weak. They are often associated with the presence of opposite magnetic 
polarities within the spatial resolution clement and a magnetic vector that is 
almost perpendicular to the line of sight, situations which naturally give rise to 
small V [2E1I2II- 

The observed Stokes U asymmetry is on average smaller than the V asym- 
metry. This is true in particular for extreme asymmertic values, i.e. 
(C / a,obs/C / s,obs)max < (V r s ,obs/^a,obs)max- Since this relation also holds at fx = 1 
it is valid for source profiles as well. Thus, Sanchez Almeida & Lites (1992) 
j2E] point out that Stokes U retains t/ a! obs/^s obs ^ 1 throughout a sunspot, 
although V s /V a > 1 is invariably achieved at the neutral line. The reason for 
the smaller maximum asymmetry lies in the fact that Stokes U senses the trans- 
verse magnetic field. Since velocities in the solar atmosphere are directed mainly 
along the field lines they generally have a small line-of-sight component when U 
has a significant amplitude. Sizable line-of-sight velocities are needed, however, 
to produce a significant asymmetry [21]. Another reason for the smaller max- 
imum U asymmetry is that, unlike Stokes V, it does not distinguish between 
oppositely directed magnetic fields. 

Thus it is not surprising that in the following analysis Stokes U provides 
tighter limits than Stokes V. Another reason is that due to the on average 
stronger observed V profiles, asymmetries introduced in U (through gravita- 
tionally introduced cross-talk from V) are larger than the other way round. 
However, we also analyse Stokes V as a consistency check. 

In order to seperate the asymmetry produced by solar effects from that 
introduced by gravitational birefringence, one strategy to follow is to consider 
large amplitude Stokes profiles only. Another is to analyse data spanning a large 
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range of fj, values, since A$ exhibits a definite fi dependence. Finally, the larger 
the number of analysed line profiles, the more precise the limit that can be set 
on A$. Better statistics not only reduce the influence of noise, they are also 
needed because for a single profile gravitational birefringence can both increase 
or decrease V s /V a and U a ,/U s . The latter may become important if the source 
profiles are strongly asymmetric. Thus a small observed V s /V a or U^/Us is in 
itself no guarantee for a small gravitational birefringence. However, since almost 
all source profiles are expected to have V s /V a <C 1, U^/Us < 1, on average we 
expect gravitational birefringence to increase these ratios. 

4 Observations and data 

Two sets of data have been analysed in the present paper. They are described 
below. 

4.1 Data obtained in 1995 

Observations were carried out from 7 - 13th Nov. 1995 with the Gregory Coude 
Telescope (GCT) at the Teide Observatory on the Island of Teneriffe. For the 
polarimetry we employed the first version of the Zurich Imaging Polarimeter 
(ZIMPOL I), which employs 3 CCD cameras, one each to record Stokes I ±Q. : 
I ±U and I ±V simultaneously . 

The recorded wavelength range contains four prominent spectral lines, Fe 
I 5247.06A, Cr I 5247.56A, Fe I 5250.22A & Fe I 5250.65A. Three of these 
spectral lines are among those with the largest Stokes amplitudes in the whole 
solar spectrum and are also unblended by other spectral lines Blending 
poses a potentially serious problem since it can affect the blue-red asymmetry 
of the Stokes profiles. By analysing more than one such line it is possible to 
reduce the influence of hidden blends and noise. Nowhere else in the visible 
spectrum are similar lines located sufficiently close in wavelength that they can 
be recorded simultaneously on a single detector. Also, compared to other lines 
with large Stokes amplitudes the chosen set lies at a short wavelength. This 
enhances A$ since it scales with 1/A. The sum of the above properties make 
the chosen range uniquely suited for our purpose. 

In order to image all 4 spectral lines of interest onto a single CCD we in- 
troduced reduction optics between the image plane of the spectrograph and the 
detectors. They produce an image-scale reduction by a factor of 3.2. The final 
spectral resolving power A/AA corresponded to 210'000. The spatial scale cor- 
responding to a pixel was 1.13" (or 860 km on the Sun). However, the effective 
spatial resolution of the data is limited by turbulence in the Earth's atmosphere, 
so-called 'seeing'. This varied somewhat in the course of the observing run, so 
that the estimated angular resolution of the observations lies between 2.2" and 
3". 

The modulator package, composed of 2 photoelastic modulators oscillating 
at frequencies of around 41 kHz & 42 kHz followed by a Clan linear polarizer, 
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was placed ahead of the entrance slit to the spectrograph, but was nevertheless 
(unavoidably) located after 2 oblique reflections in the telescope. Oblique reflec- 
tions produce cross-talk between Stokes parameters, i.e. they partially convert 
one form of polarization into another. Since we are trying to observe, or at 
least set limits on "cross-talk" between Stokes U and V due to gravitational 
birefringence we took some trouble to reduce the instrumental cross-talk to the 
extent possible. A first step was the choice of the telescope. With only two 
oblique reflections, whose relative angles change only slowly in the course of 
a year, the GCT is relatively benign compared to most other large solar tele- 
scopes. Secondly, a half-wave plate was introduced between the two oblique 
reflections. Sanchez Almeida et al. |301 13*1] have pointed out that a half-wave 
plate at that location should, under ideal circumstances, completely eliminate 
all instrumental cross-talk. To test the efficiency of the half-wave plate in sup- 
pressing instrumental cross-talk between Stokes Q, U and V we first carried out 
a series of observations of a sunspot near the centre of the solar disc both with 
and without a half-wave plate introduced in the light path. Such tests were 
necessary since the half-wave plate available at the GCT is not optimized for 
the observed wavelength. Note that at solar disc centre (/i = 1) A$ = 0, so 
that we test for instrumental cross-talk only. The half-wave plate was indeed 
found to significantly reduce instrumental cross-talk. Remaining cross-talk was 
removed during data reduction using a numerical model of the telescope that 
includes an imperfect half-wave plate (adapted from a model kindly provided 
by V. Martinez Pillet). The parameters of the model were adjusted slightly 
using the observations of a sunspot umbra located close to fi = 1. We estimate 
that the residual cross-talk after this procedure is at the level of a few percent. 
Since Stokes V, Q, U profiles generally have amplitudes of 0.1/ c or less, the 
influence of the cross-talk is of the same order as the noise, which is roughly 
1 — 2 x 10~ 3 / c , where I c is the continuum intensity. Photon noise is by far 
the largest contributor to this noise level. At this level instrumental cross-talk 
ceases to be of concern for our analysis. 

The Zimpol polarimctcrs arc unique in that they combine CCD detectors 
with a very high modulation frequency (> 40 kHz) and hence preclude distortion 
of the Stokes profiles and cross-talk between them due to seeing fluctuations. 
The other advantage of the ZIMPOL working principle is the fact that the 
fractional polarization is absolutely free from gain-table or flat field noise |28| . 
This again improves the accuracy of the profile shapes and hence the reliability 
of our results. 

A total of 106 recordings were made at different locations on the solar disc 
in an attempt to cover a large range of fi homogenously. Particular emphasis 
was placed on observations close to the limb since gravitational birefringence is 
expected to be largest for such rays. 

Since only a single sunspot was present on the solar disc during the observing 
run most recordings refer to faculae and network features, i.e. magnetic features 
with lower Stokes Q, U, V signals. In Fig.[2]a sample Stokes J, Q, U, V spectrum 
of a facular region near the solar limb is plotted. The 4 analyzed spectral lines 
are identified. These data were fully reduced following the tedious, but well- 
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Figure 2: Sample spectrum of Stokes /, Q, U and V. 
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tested procedures described by Bernasconi (1997) |3*2). 
4.2 Data set of March 2000 

In order to improve the statistics and the fj, coverage a second observing run 
was carried out in March 2000 with the Gregory - Coude Telescope of IRSOL 
(Istituto Ricerche Solari Locarno) in Locarno, Switzerland. This telescope is 
almost identical to the GCT on Tenerife and the parameters such as spectral 
resolution, noise level etc. are very similar to those of the 1995 observations. 

The next generation, ZIMPOL II polarimeter [JJHJ [5Ij was employed for the 
polarization analysis and data recording. It simultaneously records three of the 
four Stokes parameters, either Stokes /, Q, V or /, U, V on a single CCD detector 
chip. Observations in these two modes were interlaced, such that alternate 
exposures record Stokes /, Q, V and /, U, V, respectively. Exposures of the same 
Stokes parameters were then added together to reduce noise. Thus the final data 
set consists of all four Stokes parameters. The only differences with respect to 
the recordings made in 1995 are that the number of spatial pixels is reduced 
and that the noise level of Stokes V in the newer data is a factor of V2 lower 
than of Stokes Q and U, whereas Stokes Q, U and V had the same noise level 
in the earlier recordings. Due to the superior modulation scheme implemented 
in ZIMPOL II Stokes Q and U achieve a noise level of 10 _3 / c after roughly the 
same exposure time during the observations made in 2000 as during the earlier 
campaign. 

These observations were carried out on the day of the equinox, at which 
time the two mirrors producing oblique reflections of the beam ahead of the 
modulator package are oriented such that their polarization cross-talk cancels 
out. Hence for these observations the instrumental cross-talk is essentially zero 
and no further treatment of the data for this effect is required. 

The Sun was very active at the time of these observations with many active 
regions harboring sunspots and faculae present on the solar disc. Since active 
regions generally give larger amplitude Stokes signals we concentrated on ob- 
serving them. A total of 7 exposures were made. The typical seeing during 
these observations was estimated to be 2 - 3", while the spatial pixel size was 
1.13". 

5 Data analysis 

Each exposure gives us the profiles of the 4 analyzed spectral lines in Stokes 
/, Q, U and V at a set of 94 (128 in the 1995 data) positions on the solar disc. 
Once the reduction and calibration procedure is completed we select from a 
given frame those spectra for further analysis for which the S/N ratio for either 
Stokes U or V in at least one of the four spectral lines is above 12 in the 1995 
data and above 15 in the 2000 data. 

This criterion gave us a total of 4480 profiles for further analysis. A non- 
orthogonal wavelet-packets smoothing scheme was employed to enhance the S/N 
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Figure 3: Simultaneously measured profiles of Stokes V and Stokes U of the Fe 
I line at 5250.22A. 



ratio by a factor of 1.5 - 2 without significantly affecting the profile shapes |35|. 

Then the signed amplitudes of the blue and red wings, ab and a r (i.e. of 
the blue and red Zeeman cr-component) of all Stokes profiles of all 4 spectral 
lines was determined. In Fig. we plot observed Stokes U and V profiles of 
the Fe I line at 5250. 22A. r (V) and % r (L/) are indicated in the figure. Using 
these we can form the symmetric and antisymmetric parts of the Stokes V and 
U profile amplitudes, 14 = (a D + a r )/2, V a = (a D — a r )/2, U s = (% + «r)/2 and 

= («b — fl r)/2, respectively, which enter Eqs. ffify and l|27|l. In the following 
all V and U values and parameters are normalized to the continuum intensity, 
even when not explicitely mentioned. 



5.1 Profile difference analysis 

We now apply the technique outlined in Sec. III. A to our data. Due to the limited 
number of profiles and their irregular distribution over ijl (see Sect. V. B) any 
limit on gravitational birefringence will be less tight than what is achievable 
with ideal data presented in Sect. III. A. In Fig. 0] we plot 

vs. k, (28) 



< IVobsl > + < l^obsl > 



for different initial phase differences A$ between the orthogonal modes of line Fe 
I 5250. 65A. The averaging has been done over the /i values at which observations 
are available. This line is chosen, since it gives the tightest limits on k. The 
thick horizontal line represents the value obtained from observations. Obviously 
above k 2 = (13.8 km) 2 the curve obtained from theory always lies below the 
observed value, hence ruling out such k values. 
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Figure 4: The observable (Stokes V, Stokes U) mixture for initial phase differ- 
ences of Oir, 0.17T, 0.257T and 0.47T plotted vs. MAG coupling constant k. The 
horizontal solid line represents the value obtained from observations. Note that 
A$ values bigger than 0.57T give cyclic results. 



For comparison with literature values it is useful to set limits on the £q- 
parameter in Moffats NGT, since earlier work has concentrated on constraining 
this theory. Using the profile difference technique we obtain £q < (178 km) 2 
measured in the line Fe I 5250. 65A, compared with the previous tightest upper 
limit of (305 km) 2 . 

5.2 Stokes asymmetry technique 

A measure of the asymmetry of a Stokes profile is given by the ratio 5V = Vs/V^ 
respectively 5U = U a /U s [23 E]) In Fig. we plot these quantities vs. \V a \ 
and |Z7 S |, respectively. Each point in these plots refers to a Stokes profile of 
the Fe I 5250.65A line. Although our observations cover a range of \x values, 
Fig. [S] is similar to corresponding figures based on data obtained near ji = 1 
|23II24| . where the influence of gravitational birefringence vanishes for symmetry 
reasons. For large amplitudes (V a , U s ) the relative asymmetry (V s /V a , U a ,/U s ) 
is small, while for weaker profiles it shows an increasingly large spread. For the 
weaker profiles this spread is mainly due to noise as can be judged from the 
solid and dashed curves in Fig.0 which outline the la and 3a spread expected 
due to photon noise, respectively. The curves reveal that Stokes profiles with 
amplitudes (V a , U s ) smaller than one percent of the continuum intensity are so 
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Figure 5: Measured amplitude asymmetries for Stokes U and Stokes V. 



strongly affected by noise that they are of little use for the present purpose. This 
leaves us with 1966 individual profiles for further analysis. In Fig. H3 we plot a 
histogram of the number of these profiles as a function of \x. The distribution is 
uneven, beeing determined by the position on the solar disc of magnetic features 
at the times of the observations. 

The further analysis is made more complicated by the fact that a V s and a 
U a signal can be produced not just by gravitational birefringence, but also by 
radiative transfer processes acting in the dynamic solar atmosphere, as described 
in Sect. 3. To circumvent this problem we consider all profiles satisfying the 
criterion that \V a \ or \U S \ > 0.01. For all these profiles |(5V^bsl < 0.7 an d 
|<5£/ b s | < 0.6. A similar picture is also obtained at fi = 1. Hence one way to 
limit k 2 is to require |5T4 rc | < 1 and |<5f7 src | < 1 for all [i. This condition is 
strengthened by the fact that |(5V^,bsl an d l^obsl decrease with decreasing /i, 
|2ZI 03] , whereas gravitational birefringence increases towards the limb, so that 
one would expect exactly the opposite behaviour if gravitational birefringence 
had a significant effect on the Stokes V or U profiles. 

In Fig. \7\ we plot the maximum (C/ a ,src/C4,src) value predicted for each of 
the 4 spectral lines, based on all analysed data, vs. k 2 . The horizontal line 
represents log([/ a ,src/E4,src = 1), a limit above which this ratio is not observed at 
Clearly, as k 2 increases (f7 ajSrc /[/ s ,src)max initially remains almost equal 
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Figure 6: Histogram of the number of recorded profiles as a function of /i. 

to (f/ a dbs/U s obs)maxj but begins to increase for fc 2 >(lkm) 2 , becoming > 10 
at k 2 < (2 km) 2 and finally oscillating around (f/ a ,src/£/s,src)max of 100 - 1000. 
All 4 spectral lines exhibit a similar behaviour, implying that the influence of 
noise is very small. The largest effect of gravitational birefringence is exhibited 
by the two most strongly Zeeman split lines, Cr I 5247.56A, Fe I 5250.22 A, 
which also produce the largest Stokes V and U signals, while the lines with the 
smallest splitting, Fe I 5250.65 A provides the weakest limit. 

It is in principle sufficient to limit k by requiring that none of the observed 
spectral lines has (£7 a ,src/E^s,src) > 1 within the range of allowed k values. 
This gives k 2 < (0.9 km) 2 . A limit obtained similarly from Stokes V is both 
larger k 2 < (3.42 km) 2 and less reliable, since we cannot completely rule out 
Vs,src/Va,src > 1 to be present, although we expect such profiles to be very rare, 
among large amplitude profiles. 

For Moffat's NGT this technique gives l\ < (57, 1km) 2 , implying a 28-fold 
reduction in £q compared to previous work. 

An alternative test is to determine the fraction of profiles with (£/a src /[/ s src ) > 
1 (Fig. El . Fig. [S] reveals that initially no U profile satisfies the criterion, above 
k% = (1.33 km) 2 1.7 % of the profiles does. This number keeps increasing with 
k, before finally oscillating around 70% at large k. Thus 10% of all data points 
have (£/ a ,src/E/s,src) > 1 for k 2 = (1.955km) 2 , 20% for k 2 = (2.43km) 2 . We are 
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Figure 7: Maximum of (t/ a)Src /i7 S)Src ), on a logarithmic scale, values for all 4 
spectral lines vs. k. 

not aware of any solar observations of Stokes U with SU > 1 for which instru- 
mental cross-talk is negligible |2fill38| . k 2 < (2.5 km) 2 is thus a very conservative 
upper limit, deduced from this criterion. 

For £ 2 , we found that 10% of all data points have (J7 a ,src/£4,src) > 1 f° r 
£% = (74.0km) 2 and 20% for £% = (79.4km) 2 . 

6 Discussion and Conclusions 

New couplings between electromagnetism and gravity introduced within the 
framework of Metric-affine Gravity (MAG) theories lead to observable effects, 
in particular to space-time becoming birefringent in the presence of a gravita- 
tional field. By constraining the level of this birefringence we can constrain 
the strength of the coupling between electromagnetism and gravity within this 
particular framework. A birefringence of space-time is also expected in a more 
general context, so that such a constraint is also of wider significance. 

Using two techniques (the Stokes asymmetry technique and the new profile 
difference technique applied to solar data) we have imposed stringent constraints 
on the coupling constant k introduced by MAG. In order to judge how strin- 
gent the constraints imposed by these techniques are, we have also redone the 
analysis for Moffat's NGT. The new techniques improved previous limits on i 2 
given by Solanki & Haugan (1996) 4 by nearly one order of magnitude. It will 
be difficult to set much tighter limits on gravitational birefringence than those 
found here using solar data in the visible spectral range. To obtain a significant 
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Figure 8: Fraction of profiles with {U a>BTC /U BjSIC ) > 1 for NGT (top) and MAG 
(bottom). 



improvement one would need to observe at shorter wavelengths. The spectral 
line at the shortest wavelength that is strong enough to provide a hope of de- 
tecting Stokes U and V at sufficient S/N is Lya at 1216A. The maximum gain 
that one could expect relative to the current analysis is a factor of 

iVisible = 5250 = 

A L ya 1216 V ' 

Another possibility is to consider astronomical objects with a stronger gravita- 
tional field and well defined sources of polarization. In a parallel paper |39| we 
employ magnetic white dwarfs for this purpose. 
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